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Abstract 

We study the one-dimensional contact process in its quantum version using a recently pro- 
posed real space renormalisation technique for stochastic many-particle systems. Exploiting the 
duality and other properties of the model, we can apply the method for cells with up to 37 sites. 
After suitable extrapolation, we obtain exponent estimates which are comparable in accuracy 
with the best known in the literature. 
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1 Introduction 



Phase transitions out of an absorbing state form an important class of non-equilibrium critical 
phenomena fij. Models having such a transition have appeared in various areas such as surface 
chemistry pj, population dynamics ...Very recently it was even shown that the so called 
'self-organised criticality' appearing in a number of sandpile models can be related to 'ordinary' 
criticality in a class of models with an infinite number of absorbing states and a conservation law 

Hi- 

A hot topic in current non equilibrium statistical mechanics is to understand what are the possible 
universality classes that can exist for such models, and more importantly, what precisely determines 
these universality classes. Research on these questions has mainly focused on the one dimensional 
case. By now, it is clear that two main non-trivial universality classes exist. The first one is that 
of directed percolation (DP) which contains models such as the contact process, the Ziff-Gulari- 

Barshad Q model of catalysis, branching and annihilating walks with an odd offspring ||, This 

class is very robust in agreement with the DP conjecture jjj which states that all phase transitions 
out of an absorbing state in models with a scalar order parameter, short range interactions and 
no conservation laws belong to the DP universality class. Over the last decade, the existence of a 
second universality class has been clearly established. This class, known as the parity conserving 
(PC) class contains, among others, such models as the branching and annihilating walks with an 
even offspring, the monomer-dimer model §|, and a certain type of generalised contact process 
|J. The precise conditions which determine this class are however still unclear. Some authors 
argue that it is a conservation law [|l0| (the conservation of particle number modulo 2) that is the 
important factor, whereas others have claimed that it is the existence of two equivalent absorbing 
states ||. 

Another issue which is currently debated is the importance of exclusion in these models. Indeed, 
it has been recently argued that adding exclusion to a model of branching and annihilating walks 
with N >2 species of particles changes some of the critical exponents uW. 



In the light of these questions, the development of precise approximate techniques is crucial. Most 
of the current understanding of these models has come from two approaches: extensive numerical 
simulations and field theoretic renormalisation group (RG) techniques (l^, [l3|]. Both methods have 
their strong and weak points. Simulations allow the study of quite big system sizes (especially in 
d = 1). Near the critical point however, relaxation times may be quite large and one can never 
be very sure that the asymptotic time regime has been reached. Field theoretic techniques are 
very powerful but have their own difficulties. In the case of the branching and annihilating random 
walks with an even offspring, there exist two upper critical dimensions which make reliable exponent 



estimates in d = 1 very difficult [10] 



For these reasons, in the last years attention has been given to alternative approaches which can 
be called real space renormalisation methods. These techniques use the by now well known formal 
equivalence between a stochastic system and a quantum mechanical model evolving in imaginary 



time [14]. In this way, one can associate with the generator of a Markov chain a quantum Hamil- 
tonian, which can then be studied using various techniques that were originally introduced in the 
study of quantum spin chains, fermion models,. . .For some models this approach can lead to an 



2 



exact solution. As an example, we mention the relation between the asymmetric exclusion process 
and the XXZ-chain. Unfortunately, no models with a non-trivial bulk phase transition can be 
solved in this way. Yet, within the same spirit one can then use approximate techniques originally 
developped for quantum systems to study stochastic systems. The most famous of these approaches 
is certainly the density matrix renormalisation group (DMRG) ]i~5| , 16]. This technique has by now 
been adapted to stochastic systems [17, [l8|, 19 1. The method is asymptotic in time, but at this 



moment can treat only systems consisting of approximately 50-100 sites. The name DMRG is a 
bit of a misnomer since no renormalisation group flows are calculated. This may make it hard in 
some cases to get clear results on issues of universality. 



Another approach working within the same spirit was recently introduced by the present authors 
1 20]. In our real space renormalisation group technique we apply the so called standard renormali- 
sation method (SRG), also known as SLAC-approach [21], to the quantum Hamiltonian associated 
with the stochastic model. The SLAC approach was introduced in the late seventies in the study of 
lattice gauge theories and was subsequently applied to a great number of quantum spin and fermion 
systems [22]. In a previous paper we adapted this technique to the study of stochastic systems 



[20], and applied it to some exactly solvable cases. Surprisingly, in many cases exact results were 
recovered. We also applied our technique to the contact process, using a small cell of only 3 sites. 
The numerical estimates of critical properties that we obtained for this process were within 10 
percent of the best known values. In the present paper we extend our calculations for the contact 
process to as large cells as possible. Using underlying properties of the contact process we are able 
to get results for cells with up to 37 sites. To the best of our knowledge, this is a 'world record' 
for the SRG approach to quantum systems. Combined with good extrapolation techniques, these 
results allow us to determine very accurate exponents for the contact process. In fact, the results 
we obtain are of the same accuracy as those obtained from series expansions and simulations, and 
are more accurate then those coming from the DMRG. 



This paper is organised as follows. In the next section, we introduce the contact process and 
its quantum description. We also translate the duality of the model into a quantum mechanical 
language. In section 3, we give a brief outline of our real space renormalisation method. In section 4 
we show how for the case of the contact process, the RG flow can be calculated from the knowledge 
of only two matrix elements. In section 5 we present the results of our calculations for different 
system sizes, describe the extrapolation procedures and compare our results with those in the 
literature. Finally, we conclude with a discussion in section 6. 



2 Quantum desciption of the contact process 



The contact process was originally introduced as a simple model for the spreading of an epidemic 
[||]. On each site i (i : 1, . . . , N) of a lattice A there is a variable rji which can take on two values, 
referred to as A and 0. In the contact process, particles A (vacancies 0) are interpreted as sick 
(healthy) individuals. The dynamics of the model is given by a continuous time Markov chain on 
the set of all microscopic configurations rj = {r/i, . . . , ryjv}- The following processes are allowed: a 
sick person can cure (A — > 0) with rate 1 and a healthy individual can be become contaminated 
with a rate z\/2 where z is the number of sick neighbours. The conditional probability P(rj, t; t]q,0) 
that the system is in configuration r/ at time t given that it was in t]q at time t = then obeys the 
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master equation 



j = -}^H{r],7])P{r] ,t;r] ,0) (1) 

n' 

where the 2^ x 2 N matrix H, the generator of the Markov chain, depends on the transition rates, 
i.e. on A. 

The master equation (|]) is formally equivalent to a Schrodinger equation in imaginary time. It has 
therefore become common to introduce a quantum mechanical notation for a stochastic process. 
This mapping of a stochastic system on to a quantum mechanical one is by now quite standard and 
we will not discuss it here in any detail. We only give a brief review that fixes also the notation 
that will be used further. To each configuration r\ a state vector \rf) = ®iL±\r]i) is associated. The 
vectors |r/j) form the basis vectors of a two-dimensional vector space. It is then natural to use a 
spin 1/2 language to describe this vector space. As usual, a particle (vacancy) will be associated 
with spin down (up) p3| . 



Next, we also associate a vector \P(t)) with the conditional probabilities P(rj,t;r]Q,0) such that 

\P(t))=J2 P ( r l,t;r ] o,V)\v) (2) 
v 

Using this notation, the master equation (Q) is rewritten as 

d\P(t)} 



dt 



H\P(t)) (3) 



From now on, we will refer to the matrix H as the Hamiltonian of the stochastic system. For 
processes with transition rates that involve only nearest neighbours sites (such as is the case for the 
contact process), H can be written as a sum of local hamiltonians /ij,j+i that act only on nearest 
neighbour sites 

H = li ® . . . ® <g> h iii+ x eg) U+2 <S>... In (4) 

i 

In the particular case of the contact process, we have 

h i)i+1 = (raj - sf) (g) + ^ (vi - s~) (g> n i+ i + rii® {v i+l - s~ +1 ) (5) 
where the matrices v,n,s + and s~ are given by 



The formal solution of the master equation (^) is 

\P(t)) = e- m \P(0)) 

Because of the properties of stochastic matrices, there is always a zero eigenvalue and the real 
parts of the other eigenvalues of H are never negative. Therefore, asymptotically, for t — > oo 
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the behaviour of \P(t)) is determined by the properties of the ground state(s) of the quantum 
Hamiltonian. In our RG approach, we study the critical behaviour of the stationary state of the 
stochastic system by applying a real space renormalisation technique originally developped to study 
ground state properties of quantum systems. 

Expectation values of physical quantities can also easily be rewritten in terms of the quantum 
notation. With each physical quantity T (such as the density of particles, correlation functions, 
. . . ) we can associate a quantum mechanical operator F (with matrix elements (rj\F\r]') = J'ijjjSjm 1 ) 
such that the expectation value of T 

n 

can be rewritten as 

(T)(t) = (s\F\P(t)) = (s\Fe- m \P(0)) 
where we have introduced the short hand notation 

(s\=J2(v\ (7) 
v 



The contact process has a property known as duality |]24|]. This notion was first introduced in the 
probabilistic study of interacting particle systems and should not be confused with the concept of 
duality from equilibrium statistical mechanics. Before proceeding with the renormalisation group 
study of the contact process, we show how this duality can be derived in the quantum mechanical 
language. To the best of our knowledge, this formulation of the duality of the contact process 
has not yet appeared in the literature. Moreover, we will use it to simplify the renormalisation 
calculations (see section 4). 

We begin by introducing the 2x2 matrix d 




and dk = li <S> I2 • ■ ■ <8> d <2> . . . ® ljv, VA; (with d on the k-th site). It is then easy to check the 
fundamental property 

hi = (4 <g) d k+ i)h k (d k ® d k+1 y l 

where T stands for transposed. Furthermore one has d k n k d~^ 1 = v k — sf. From these relations and 
defining D = ® k d k one can obtain a useful expression for the expectation value of n k which gives 
the density of particles at site k 

(n k )(t) = {s\n k e- Ht \P(0)) 

= {s\D- 1 Dn k D- 1 Der Ht D- 1 D\P(0)) (8) 

If, |0) and \N) denote respectively the completely empty and completely full configuration one has 

(D-y\ s ) = |o> 

D\L) = |0) 
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Therefore, if we take as initial condition \P(0)) = \L), @ becomes 

<n*>(t) = (0|(^-4)e-^|0) 
= (0\e- Ht (v k -s^)\0) 

= l-(0\e- Ht \k) (9) 

Here (A;) is the state with only one particle, at site k. The matrix element (0|e~^*|A:) appearing in 
@ gives the probability that starting at t = with one particle at site k no particles are left in 
the system at time t. If we introduce the survival probability Pk{t) as the probability that if the 
system is initially in the state \k), there are still particles in the system at time t we finally get 

(n*)(t) = P k (t) (10) 

This is the duality relation of the contact process. It says that the density of particles at site k 
when starting from a completely full lattice is the same as -Pfc(i). For t — > oo, one expects that the 
steady state value of the density becomes independent of the initial condition and hence we obtain 

{n k )st = Pk,st (11) 
Here, the subscript st denotes the stationary state value. 

(In the rest of this paper we will drop the direct product symbol to shorten the notation, we will 
also drop all unity operators, their presence is always implicitly assumed) 



3 Standard renormalisation for stochastic systems 



We briefly review the use of the standard real space RG for stochastic system as introduced in a 
previous paper [20]. For more details we refer to that work. 



As usual in real space RG approaches, the lattice is divided into cells, each containing b sites. In 
the case of a one dimensional system, we can regroup the terms in the Hamiltonian H @ to write 

H = ^2(Ho,a + V a , a +l) (12) 
a 

Here a labels the cells, PiQ a contains the intracell terms of H and V a ^ a +i the intercell interactions. 
Next, i?o,« is diagonalised exactly. For simplicity, we now only consider the case where the ground 
state of -ffo.a is doubly degenerated. We then have two right and two left ground states of Hq iCI 
denoted as |si) a , \s2) a and a (si\, a {s2\ which we can normalise as a {si\sj) a = 5ij. We consider one 
of these states as representing a 'cell vacancy' state |0) Q , the other as a 'cell particle' \A) a state. 
These states are used to construct renormalised lattice configurations \rf) = which span a 

2 N / b dimensional subspace VV of the original state space. 

The renormalisation transformation is now performed by projecting the original Hamiltonian onto 
W. This is done by means of the matrices 

T l = ^2\ e v')( r l'\ T 2 = ^2\ri'){e v '\ (13) 

rf rf 
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where \e v >) are the vectors of the standard basis of W. Because of the normalisation of the ground 
states that we have choosen T1T2 = 1, the identity operator on W. Finally, T\ and T2 are used to 
calculate the renormalised Hamiltonian H' as 



H' = TiHT 2 (14) 

When the ground state of the intracell parts i?o,o is doubly degenerated (as we assumed), it is easy 
to show that H' is again stochastic. 



If we collect the rates appearing in H in a vector w, (14) defines a mapping in the parameter space 
w' = f{w). From this mapping we can determine fixed points, critical exponents, expectation 
values in the ground state, . . .[] 



To fix ideas, let us assume that the equations w' = f{w) have a non-trivial fixed point at w* , 
with one relevant scaling field (which in linear approximation is proportional to Awi = w\ — w\) 
whose scaling dimension is y wi . From standard RG theory it then follows that near criticality the 
correlation length £ will diverge as £ ~ lAu^l" 1 - with 

vx = 1/ywi (15) 



In order to determine the order parameter exponent (5 we need first to explain how the particle 
density in the stationary state c s t{w) can be calculated within our RG scheme. When the system 
is in the ground state |sj), and assuming translational invariance, this density is given by 



c s t{w) = (s\n k \si(w)) 



(16) 



where we have now explicitly indicated the dependence of the ground state on the transition rates 
w. Under the renormalisation this expectation value transforms as pOj] 



Cstiw) = a(w)c st (w') 



(17) 



Here we assumed, as will turn out to be the case for the contact process, that the renormalised 
particle operator n' k = T\n^T2 is proportional to n&, i.e. n' k = a(w)rik- The relation ( |l7|) can be 
iterated along the RG flow, and hence the density of particles can be obtained as an infinite product 
if one knows the density at the (trivial) fixed point w* which attracts w 



c s t(w) 



u=o 



(18) 



In principle, other expectation values can be calculated in a similar way. 



To conclude this section we show how (|l7|) can be used to calculate the exponent (3. Near w* we 
get for the singular part of c s t 

c st (A Wl ) = a(w*)c st (by™i Awi) (19) 

One of the rates appearing in the Hamiltonian can always be taken equal to one. This is nothing but a fixing of 
the timescale. The corresponding rate in the renormalised Hamiltonian is not necessarily one, but we then divide H' 
by this renormalised rate - the effect of this division is included in the mapping /. 



7 



We write 

a(w*) = bP/ v ± 

and get from (IS) and (|i~5|) 

c st (A Wl ) ~ (A Wl f 
which justifies (|20|). Hence, j3 can be obtained from a(w*). 



(20) 



(21) 



4 Renormalisation of the contact process 



This section is a more technical one, and shows explicitly how we calculate the RG flow for the 
contact process. 

We start by dividing the lattice in blocks of length b. Within a cell the terms of the Hamiltonian 
[H|) are regrouped. We therefore introduce the following short hand notations 



hi 



rii-s 
A 



(vi - ^ )n i+ x + ni(v i+ i - s i+1 ) 



(22) 
(23) 



which are respectively the generators of the processes A — > and A + % ^ A + A, $ + A ^ A + A. 
Notice that each of these terms itself has the property of duality. 



It is now important to remark that the regroupment in intra- and intercell parts of the Hamiltonian 
is not unique. A natural attempt is to take into i?o,a a U the terms that act on the sites inside 
the cell, hence take the intracell Hamiltonian as that of a contact process for a system of b sites 
with open boundary condition: i/o,a = Ya=i ^ai + Yl\~ l ^ai (here the first subindex labels the 
cell, while second one indicates the site in the given cell). This choice however is not suitable for 
us since in this case the intracell Hamiltonian has only one groundstate which is the trivial empty 
lattice \si) a = ®i=i|0)a,i- 



As we argued in the previous section, Hq^ should have two ground states, one representing the 
effective vacancy of cell a and one representing the effective particle. It is the second one that is 
missing. To solve this problem we 'force' i?o,a to have an active ground state by removing hj on 
the central site of the cell. This resembles the so called self dual renormalisation group introduced 
earlier in the study of quantum models such as the Ising model in a transverse field [25]. 



From now on, we choose b odd, b = 2n — 1 and take 

n-l b 6-1 

ffo,« = XXi+ E ( 24 ) 

i=l i=n+l i=l 

The operator h\ n which we removed from Hq^ is then of course added to V^, jQ ,+i concluding the 
regroupment of the terms occuring in the Hamiltonian. 
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The hardest part of the RG is the calculation of the two right and left groundstates of Ho t0l for 
b as large as possible. Analytically this can only be done for b up to 5 sites. To make a reliable 
extrapolation for b — > oo possible we need a good numerical algorithm to get to larger b. Before we 
turn to this point in the next section, there are a few analytical considerations that can be made 
and that what will turn out to be extremely useful in studying large cell sizes. 

We start with the two right ground states |si) Q , \s2) a of -Ho, a as defined in (p4[). The first one is 
again trivial |si) Q = <8>f = i|0) a ,i- To get an idea of how the active ground state \s2) a looks like we 
rewrite ( |2"4| ) as 

H^ a = H l a + H r a (25) 

where 

n—1 n— 1 

H a = X h ct,i + X h a,i 
i=l i=l 
b 6-1 

#a= E <i + E<i ( 26 ) 

i=n+l i=n 

Physically {H l a ) is the stochastic generator of the contact process on a lattice of n sites without 
the process A — > on the right (left) site. These operators have again a trivial right ground state 
(the empty lattice) and a non-trivial (active) right ground state. To get a better grasp on the latter 
one, we first notice that H r a has no (A — > 0)-term on site (a,n). Since this is the only reaction 
destroying ^4's on that site, there are no transitions possible from configurations with an A on site 
(a, n) to configurations with no particle at (a,n). Denote the subspace spanned by the former 
configurations by V. Then, H r a defines a stochastic process on V, implying that H T a must have a 
ground state in this subspace. This ground state clearly cannot be the empty lattice since this is 
not an element of V. We therefore conclude that this state is the active ground state of H r a and 
that it has with probability 1 a particle A at site n. The same can be said for H l a . As a consequence 
we can write the active ground state as 

\Wa ® \A) a , n for H r a 

\A) a , n <g> \i/j) l a for H\ (27) 

where \ip) r a and \ip) l a are states on a lattice with n — 1 sites. It then follows from (|26| ) that 
|?/>)£, (8> |^4)o, in <S> \ip) l a is the active ground state of Ho ia . Hence, if we can find the right ground state 
of the 'contact process' determined by H^ on a lattice of n sites we can construct the ground state 
of the process defined by Hq^ on a lattice of b = 2n — 1 sites (since l^)^ can easily be obtained by 
reflection once \ip) r a is known). In conclusion, we have 

|S1>Q = ®Ll|0)a,i 

|*2>« = Wa ® \A) a ,n ® ( 2 §) 

For the left groundstates of Hq^ we always have the trivial ground state a (s\ = X^ 7 /! an d a 
non-trivial one. To find the latter, we exploit the duality of Hq )0[ . Indeed, from this duality 
= BHq^B' 1 , it follows that for any right groundstate \sk) a of Ho a , (B\sk) a ) T is a left 
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ground state. Denoting £ (</>| = ) we have 



[s 2 \ = a(s\ ~a (si| (29) 



We take a (s2| of this form because of normalisation reasons. Our choice garantees that a — 

Si j, which we need to conserve stochasticity as explained in the previous section. 

To conclude the calculation we perform the renormalisation transformation (|l4|) on H. Without 
any further information on |^}^, we know that the ground states (^8||2^) of H$ a are properly 
normalised and have left-right symmetry and moreover we know the state of the central site. Using 
these three properties, it is straightforward to show that the renormalised Hamiltonian H' contains 
the same terms as H, there is therefore no proliferation of interactions, and that the RG equation 
for the rate A is of the form 

v 2 

X' = A— ^ (30) 

where 



v = (s\n at i\(ip) r a ® \A) a>n ) 

W = C a (<f>\®a,n mstm«®\AUn)= r MWa =MB T W a (31) 

This means we can generate the RG map for the contact process with cell length b = 2n — 1 by 
calculating two matrix elements in the right ground state of the contact process H r a on a lattice of 
only n sites. In this way it is possible to perform the RG for rather large cell sizes. Moreover for 
each cell size, the calculations that have to be performed are rather limited. 



5 Results 



We applied the RG procedure described in the previous sections to the contact process for block 
sizes b = 3,5,.... For each size b = 2n — 1 we had to calculate the non-trivial ground state of 
the non-symmetrical 2 n x 2 n matrix K r a . Analytically, this was only possible for 6 = 3 and 6 = 5. 
For larger block sizes we turned to numerical diagonalisation methods, in particular the Arnoldi 
algorithm [26|. This algorithm is designed to calculate eigenvalues and eigenvectors of an extreme 
part of the spectrum (in our case the low lying part) of large non-symmetrical matrices. When 
the algorithm converges, it produces very precise estimates. Since for stochastic systems we know 
the value of the ground state energy exactly, we have a reliable criterium to decide on convergence 
and hence a very powerful diagonalisation tool. Using this method, we were able to perform the 
RG-calculations up to 6 = 37. 



For each 6-value the location of the critical point A c , and the critical exponent v± were calculated 
using the methods explained in sections 3 and 4. In order to determine also the exponent ratio 
(3/v± we need to calculate the quantity a(A) at the critical point (see (^)). This requires the 
calculation of some extra matrix elements. Our results for A c and the two critical exponents are 
given in table 1. 
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We have extrapolated the results of the RG calculations using the BST-algoritm [27], which is 
known to be a good tool to extrapolate finite lattice data 28]. The results are also included in 
table 1. 



In table 1 we also compare our extrapolations with those that can be found in the literature and 
which are based on a variety of other techniques. The results in the second row were obtained from 
a numerical diagonalisation of the Hamiltonian for the contact process on finite lattice. In that case 
the exponent j3 was not calculated, but we used a scaling relation |lj to obtain this exponent from 
estimates of v± and the exponent 5. The estimates in the three last rows are not for the contact 
process itself, but for other models in the same universality class. 
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11 

13 

15 

17 

19 

21 

23 

25 

27 

29 

31 

33 

35 

37 



oc 



series expansion 



diagonalisation [31] 
simulations [3C] 



series expansions 
DMRG jjD| 



3.228740192229 
3.232841532095 
3.236622341324 
3.240001893307 
3.243002363779 
3.245669884031 
3.248051604747 
3.250189366461 
3.252118590397 
3.253868772889 
3.255464377754 
3.256925727824 
3.258269778268 
3.259510752614 
3.260660654464 



3.2982(2) 



3.29785(2) 
3.29792(2) 



1.100222670443 
1.099704726572 
1.099306840428 
1.098993499409 
1.098741258486 
1.098534370274 
1.098361971472 
1.098216363626 
1.098091954504 
1.097984591005 
1.097891127814 
1.097809141091 
1.097736734033 
1.097672401665 
1.097614935125 



1.09682 (2) 



1.0969(1) 
1.09681(1) 
1.09684(1) 
1.096854(4) 
1.08(2) 



.300770659640 
.291239313449 
.284829626291 
.280220582724 
.276745857289 
.274032449042 
.271855062319 
.270069478446 
.268579043181 
.267316523888 
.266233690895 
.265295037853 
.264473840447 
.263749596182 
.263106311367 



0.2534 (4) 



0.2520(1) 
0.256(1) 
0.25208(1) 
0.252072(11) 
0.249(3) 



Table 1: Critical parameters as calculated by the RG method for a block of b sites, together with 

the results coming from other approaches. 



As can be seen, the results of our RG technique compare very well with those of the other tech- 
niques, especially for the location of the critical point and the correlation length exponent. The 
value of (3/u± is somewhat less precise. However, in comparison with the standard of real space 
renormalisation calculations, the current results must be considered as extremely precise. 
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6 Conclusions 



In this paper we have applied a real space renormalisation group technique, originally developped for 
quantum systems, to the contact process. Using some analytical properties, such as the duality of 
the process, we have been able to carry out the renormalisation for rather large cell sizes. Together 
with a suitable extrapolation, this has yielded estimates for critical properties that are of very high 
accuracy. 



In our previous paper [20] we applied the technique to simple reaction-diffusion processes which 
don't have a phase transition. There we showed that in some cases our RG technique was able to 
reproduce exact results. Combining the results of the two papers, we believe that it is fair to say 
that the technique is able to give accurate results for the stationary state of stochastic systems with 
one type of particle (or stated otherwise, in which the variable at each site can be in two different 
states), when only nearest neighbour interactions are involved. Of course, it may be so that for any 
particular model some 'cooking' is necessary in order to obtain good results. But that is a quite 
general limitation of real space renormalisation approaches. 



At this moment, there are two obvious directions in which to develop this RG method. Firstly, one 
may consider processes in which more then two particles are involved. One can think for example 
of the process + A + — > A + A + A that appears in the branching and annihilating walks with an 
even offspring, a model that belongs to the PC university class. In that case the Hamiltonian of the 
process contains three site interactions. As argued in our previous paper, it is necessary to extend 
the current RG procedure to higher order, to be able to obtain a renormalised Hamiltonian with 
three site interactions (in [ 33 1 such a higher order extension of the standard RG method is discussed 
in the context of quantum spin chains). This higher order extension requires the knowledge of all 
the eigenvalues and eigenvectors of the cell Hamiltonian, which severely restricts the cell sizes 
that can be studied, since algorithms such as the Arnoldi or Lanczos procedures are only able to 
give good estimates of low lying eigenvectors and eigenvalues. Moreover, in stochastic systems, 
eigenvalues can be complex, which in turn can give rise to parameter flows which are complex. We 
have tried this kind of higher order technique in a preliminary study of a non-equilibrium Ising 



model [34] whose transition is believed to be in the PC class, and which also has a duality [p5| . 
In that calculation we encountered this problem of complex eigenvalues and it is at this moment 
unclear to us how to proceed in this direction. 



A more promising approach extends the techniques introduced here to models with several types 
of particles, or with more then two states per site. If one restricts again the interactions to be of 
nearest neighbour type there are no fundamental problems to apply our RG technique. Several of 
the interesting processes mentionned in the introduction belong to this class of models. One may 
think of branching and annihilating walks with two types of particles and exclusion, or the model 
originally studied by Van Wijland et al. p6| , . . . 

Another very interesting model of this latter type was recently introduced by Hinrichsen ||. It is 
a generalisation of the contact process in which at each site there can be n 'empty' or non-active 
states. For n = 1, the model coincides with the contact process studied here. For n = 2, the 
generalised contact process is believed to be in the PC universality class, whereas for n > 2 the 
critical behaviour has not been determined yet Q . Because this model is a natural extension of the 
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contact process, and given the success of our RG method for that model, we believe it is an example 
of an interesting model that could be studied succesfully with our approach. Unfortunately, the 
model has no obvious duality. Moreover, since at each site, the system can be in n + 1 states, the 
calculations will by necessity be restricted to smaller b values. Nevertheless, for n = 2 it should 
still be possible to reach cell sizes b ~ 20 — 25. In this way, we hope it will be possible to obtain 
rather accurate exponent estimates for the PC universality class. We plan to present results of an 
RG study of the n = 2 Hinrichsen model in a forthcoming paper. 
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